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Abstract 

We write exact renormalizat ion-group recursion relations for a ferromag- 
netic Ising model on the diamond hierarchical lattice with an aperiodic distri- 
bution of exchange interactions according to a class of generalized two-letter 
Fibonacci sequences. For small geometric fluctuations, the critical behav- 
ior is unchanged with respect to the uniform case. For large fluctuations, 
the uniform fixed point in the parameter space becomes fully unstable. We 
analyze some limiting cases, and propose a heuristic criterion to check the 
relevance of the fluctuations. 



I. Introduction 



The experimental discovery of quasi-crystals motivated many investiga- 
tions of the effects of geometric fluctuations produced by different types of 
aperiodic structures. There are several specific results for phonon and elec- 
tronic spectra of linear problems.0 There are also some results for more 
difficult non-linear problems, as the analysis of the critical properties of a 
ferromagnetic Ising model on a layered lattice with aperiodic exchange inter- 
actions along an axial direction.! 

In a recent publication,! Luck presented a detailed study of the critical 
behavior in the ground state of the quantum Ising chain in a transverse field 
(which is known to be related to the transition at finite temperatures of the 
two-dimensional Ising model). The nearest-neighbor ferromagnetic exchange 
interactions are chosen according to some (generalized) Fibonacci sequences. 
The geometric fluctuations are gauged by a wandering exponent u associ- 
ated with the eigenvalues of the substitution matrix of each sequence. The 
critical behavior remains unchanged (that is, of Onsager type) for bounded 
fluctuations (small values of lu). Large fluctuations induce much weaker sin- 
gularities, similar to the case of a disordered Ising ferromagnet. 

In this paper, we take advantage of the simplifications introduced by 
a hierarchical lattice to give another example of the effects of geometrical 
fluctuations on the critical behavior of a ferromagnetic model. It should be 
pointed out that hierarchical lattices have been widely used as a toy model 
to test critical properties. In some cases the (approximate) Migdal-Kadanoff 
renormalization-group transformations on a Bravais lattice are identical to 
the (exact) transformations on a suitable hierarchical structure.0 There are 
rather detailed studies of these exact transformations for the Ising model on a 
variety of hierarchical structures.! There are also some recent investigations 
of an Ising spin-glass on the diamond hierarchical lattice.! 

This paper is organized as follows. In Section 2, we make some comments 
on Fibonacci sequences and introduce a nearest-neighbor ferromagnetic Ising 
model on a diamond hierarchical lattice with q branches. To simulate a lay- 
ered system, the exchange interactions {J a > and Jb > 0) are distributed 
along the branches according to the same aperiodic rule. For a large class of 
two-letter Fibonacci sequences, we write exact renormalization-group recur- 
sion relations in terms of two parameters, Xa = tanhi^, and Xb = tannics, 
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where Ka = I3Ja, and Kb = (3Jb, and (3 is the inverse of the temperature. 
In Sections 3 and 4, we analyze two typical and distinct examples. We ob- 
tain the fixed points and the flows of the recursion relations for these specific 
cases. In the first example, for small to, the critical behavior is unchanged 
with respect to the uniform case {J a = Jb)', the physical fixed point in pa- 
rameter space is a saddle point. In the example of Section 4, however, the 
fluctuations turn the physical fixed point fully unstable. We also analyze 
some limiting cases. For q = 1 (Ising chain), there is no transition at finite 
temperatures. The fixed point at zero temperature, however, may change its 
character depending on the value of the wandering exponent. In a particular 
infinite-branching limit (q — > oo), we obtain rather simple results. In Section 
5, we present some conclusions, as well as a heuristic adaptation for the di- 
amond hierarchical lattice of Luck's criterionl for the relevance of geometric 
fluctuations. 



II. Definition of the model 



Consider a particular two-letter generalized Fibonacci sequence given by 
the substitutions 

A^AB, 

B -> A A. 1 1 

If we start with letter A, the successive application of this inflation rule 
produces the sequences 

A -f AB ->• ABAA -> ABAAABAB (2) 

At each stage of this construction, the numbers and N B , of letters A and 
B, can be obtained from the recursion relations 

with the substitution matrix 

m = (1S). (4) 

The eigenvalues of this matrix, Ai = 2 and A2 = — 1, govern most of the 
geometrical properties. In a more general case (that is, for a more general 



3 



rule), at a large order n of the construction, the total number of letters is 
given by the asymptotic expression 



N n = N n + N n „ A « (5) 

where Ai > | A2 1 - The smaller eigenvalue governs the fluctuations with respect 
to these asymptotic values, 

AN n ~ AA^ = AiV™ ~ |A 2 | n . (6) 

From these equations, we can write the asymptotic expression 

AN ~ N", (7) 

with the wandering exponent 

U = W (8) 

Now we consider a nearest-neighbor Ising model, given by the Hamilto- 
nian 

W = - J ijWj> ( 9 ) 
(ij) 

with the spin variables crj = ±1 on the sites of a diamond hierarchical struc- 
ture. In Fig. 1, which is suitable for the period-doubling Fibonacci rule of 
Eq. ([!]), we draw the first stages of the construction of a diamond lattice with 
a basic polygon of four bonds (ramification q = 2). As indicated in this 
figure, we simulate a layered system by the introduction of the interactions 
J a > and Jb > along the branches of the structure. Considering the 
elementary transformations of Fig.2, and using the rules of Eq. ([!]), it is 
straightforward to establish the recursion relations 



tanh K' A = tanh 



2tanh _1 (tanhiT A tanhir B )| , (10) 



and 



tanh K' B = tanh 



2 tann' 1 (tanh 2 ^)! , (11) 



where Ka,b — PJa,b- in this particular case, these equations can also be 
written as 

2xa%b 

1 -|- X j^X q 



*'a = T^S!^ (12) 
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and 

X'r, = 



*'b = ^6r, (13) 



where 

xa,b = tannics- (14) 

In Section 3, we analyze the fixed points associated with these simple recur- 
sion relations. 

We can use similar procedures to consider much more general Fibonacci 
rules. To avoid any changes in the geometry of the hierarchical lattice, in 
this paper we restrict the analysis to period-multiplying Fibonacci inflation 
rules of two letters (the largest eigenvalue of the substitution matrix, Ai, 
gives the multiplication factor of the period). For example, let us consider 
the substitutions 

B -> A k+ \ (15) 
where k, I > 1 are two integers. From the substitution matrix, 

M=(f* + '), (16) 

we have the eigenvalues, Ai = k+l and A2 = —I, and the wandering exponent, 
u = In 1 1 In {k + 1) . Now we consider a diamond lattice with q branches and 
a basic polygon of q (k + I) bonds. The exchange interactions are chosen 
according to the general Fibonacci rule of Eq. flT5|). We then write the 
recursion relations 



tanhi^ = tanh gtanh 1 (tanh^Atanh^s) , (17) 

and 

tanh K' B = tanh [q tanh _1 (tanh k+l K^j ] . (18) 



III. Irrelevant fluctuations 



Consider again the Fibonacci inflation rules given by Eq. (fj). From the 
eigenvalues of the substitution matrix, Ai = 2 and A2 = — 1, we have uj = 0. 
For the branching number q = 2, the recursion relations are given by Eqs. 
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( |T2"D and (|13|). The fixed points and some orbits of the second iterates of 
this map are shown in Fig. 3. Besides the trivial fixed points, there is also a 
non-trivial fixed point, given by 

x* A = x* B = 0.543689..., (19) 

which comes from the solution of the polynomial equation 

x\ + 2x\ - Ax\ + 1 = 0. (20) 

The linearization about this uniform fixed point yields the asymptotic ex- 
pression 

where 

M T =fi I) (22) 



is the transpose of the substitution matrix, and 

C = 1 ~ ^ Xa) = 0.839286.... (23) 
2x* A 1 ' 

The diagonalization of this linear form gives the eigenvalues 

Ax = CAi = 2C= 1.678573..., (24) 

and 

A 2 = CX 2 = -C= -0.839286.... (25) 

As Ai > 1 and — 1 < A 2 < 0, the fixed point is a saddle point with a flipping 
approximation (in Fig. 3, we draw the trajectories of some second iterates 
of this map). Moreover, if we make J a = Jb, it is easy to see that the 
same eigenvalue Ai characterizes the (unstable) fixed point of the uniform 
model (see the diagonal flow in Fig. 3). Therefore, in this particular example 
the geometric fluctuations are unable to change the critical behavior of the 
uniform system. 

It is not difficult to check that the same sort of behavior (saddle point; 
largest eigenvalue associated with the uniform system) still holds for all finite 
values of the branching number q of the diamond structure. In the limit of 
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infinite branching (q — > oo, Ka,Kb — > 0, with g 2 AOi and g 2 K_e fixed), the 
recursion relations are particularly simple, 

Va = VaVb, (26) 

and 

Vb = vl (27) 

where = Q 2 Ka and ys = q 2 K B - The linearization about the uniform fixed 
point, y* A — y* B — 1, yields the relations 

(28) 

which correspond to a limiting case of Eq. (pi]), with C — > 1. As the limiting 
eigenvalues are given by Ai = Ai = 2, and A 2 = A2 = —1, a linear analysis 
is not enough to check the (flipping saddle-point) character of this marginal 
case. 

Another particular case of interest is the simple Ising chain (q = 1). The 
recursion relations are given by 

x' A = x A x B , (29) 

and 

x b = x Ai (20) 





with the same form of Eqs. (|26[ ) and ([27D, but with the parameters xa and 
xb given by Eq. (0). As shown in Fig. 4(a), the zero-temperature fixed 
point displays the character of a saddle point. As there is no phase transition 
at finite temperatures, it cannot be reached from physically acceptable initial 
conditions. 

IV. Relevant fluctuations 

To give an example of relevant fluctuations, consider the generalized Fi- 
bonacci substitutions, 

i^f^ (3D 
B -> AAA. y ' 
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From the substitution matrix, 



M 



1 3 

2 



(32) 



we have the eigenvalues, Ai = 3 and A2 = —2, and the wandering exponent, 
u = In 2/ In 3 = 0.630929.... For a general branching number q, we can use 
Eqs. (|17D and ([Of), to write the recursion relations 



and 



tanhi^ = tanh gtanh 1 (tanhi^Atanh 2 .?^ 



tanhi^ = tanh gtanh 1 (ta,nh 3 KA 

53) and 



(33) 
(34) 



For the particular case q — 2, Eqs. ([53 
relations 



x . 



and 



x 



B 



1 ~\~ Xa%% 



2xj 
1 + x\ ' 



) reduce to the simple 
(35) 

(36) 



where the parameters xa and xb are given by Eq. (|Hp. In Fig. 5, we show 
the fixed points and some second iterates associated with these recursion 
relations. The nontrivial fixed point, given by 



X * A = X * B = 0.786151. 



(37) 



comes from the physical solution of a polynomial equation. The linearization 
about this uniform fixed point yields the relations 



Ax' A 
Ax', 



CM 1 



B 



Ax a 
Ax B 



(38) 



where the substitution matrix is given by Eq. (|32|), and 



C 



2 (x 


* \2 
A) 


1 - (x* A f 




1 + (x A f 


2 



0.618033. 



(39) 
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From the diagonalization of this linear form we have the eigenvalues 



A 1 = CAi = 1.854101..., (40) 

and 

A 2 = C\ 2 = -1.236067.... (41) 

The absolute values of Ai and A 2 larger than unit indicate the relevance of 
the geometric fluctuations. The uniform fixed point is fully unstable (and 
there should be no transition for J a ^ Jb)- 

Again, it is not difficult to check that the uniform fixed point remains 
unstable for all values of the branching number q. In particular, for q — > 
oo, and Ka, Kb — > 0, with q x l 2 K.A and q l l 2 KB fixed, we have the limiting 
recursion relations, 

Va = VaVb, (42) 

and 

Vb = V% (43) 

where ua = q 1 ^ 2 KA and ys = q l ^ 2 KB- The linearization about the uniform 
fixed point, y* A = y* B = 1, yields the relations 

(44) 

with the substitution matrix given by Eq. (|32"D. We then have the limiting 
eigenvalues, Ai = Ai = 3, and A 2 = A 2 = —2, which confirm the unstable 
character of this fixed point. 

For the particular case of the linear chain (q = 1), we have 

x A = xax\, (45) 

and 

x' B = x\. (46) 

In agreement with the lack of a phase transition, the fixed point at zero 
temperature is unstable [see Fig. 4(b)]. 
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V. Conclusions 



From the analysis of the exact renormalization-group recursion relations 
associated with an Ising model on a variety of hierarchical diamond struc- 
tures, we show that aperiodic fluctuations of the ferromagnetic exchange 
interactions may change the character of a (uniform) fixed point. In a par- 
ticular example, with a small wandering exponent, the fluctuations are irrel- 
evant. In this case, the critical behavior is still characterized by the same 
exponents of the corresponding uniform system. In another example, how- 
ever, stronger geometric fluctuations turn the physical fixed point unstable 
in the parameter space. Even in one dimension, although there is no phase 
transition at finite temperatures, we show that the geometric fluctuations 
change the stability of the uniform fixed point at zero temperature. 

As in the work of Luck,§ it should be interesting to devise a general 
criterion to gauge the influence of the geometric fluctuations on the critical 
behavior of the Ising model on the diamond hierarchical lattice. For a large 
lattice, with q branches, the total fluctuation A J in the exchange interactions 
should be proportional to AiV. Thus we can write the asymptotic relation 

A J ~ AN ~ iV w ~ L^, (47) 

where L is a measure of the total length. The critical temperature, T c , should 
be proportional to the total value of the exchange (that is, to L q ). We can 
then define a reduced temperature, t = (T — T c ) /T c , whose fluctuations are 
given by the asymptotic form 

Tqu) 

5t~ — = L q ^- l \ (48) 

In the neighborhood of the critical point, we have L ~ £ ~ t~ u , where £ is a 
correlation length, and v is a critical exponent. Thus, we can write 



5t t~ vq( - u ~^ 
t t 

with the exponent 



(49) 



<f> = -vq(u-l)-l. (50) 

If cf) > 0, the fluctuations are irrelevant. If < 0, however, the critical 
behavior is changed drastically. To calculate a suitable value of v, we consider 
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the largest eigenvalue of the diagonal form about the physical fixed point, 
and write the usual renormalizat ion-group expression 



Ai = CAi = b yi = b l ' u . (51) 

As the largest eigenvalue Ai of the substitution matrix gives the multiplica- 
tion factor of the Fibonacci rule, we write 

b=(X 1 ) 1/q . (52) 
From Eqs. ( |5lD and (|52|) , we have 

InAi 



" glB(C7A0 
Inserting into Eq. (150), we finally have 

InAi 



(53) 



r " In (CXj 
The fluctuations are relevant for 



(u - 1) - I. (54) 



w > -LTaT' (55) 

where the prefactor C can be obtained from the linear analysis of the unstable 
fixed point of the uniform (J a = Jb) system. Although the existence of 
some counterexample should not be ruled out, the validity of this criterion 
has been confirmed by the application to a fair number of cases (including 
the examples of Sections 3 and 4). In the limit of infinite branching, the 
prefactor C tends to unit, and the criterion of relevance is reduced to the 
Pisot condition^ uj > 0. 
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Figure Captions 



Fig. 1- Some stages of the construction of a diamond hierarchical lattice 
with q = 2 branches. Letters A and B indicate the ferromagnetic exchange 
interactions (J a > and Jb > 0). 

Fig. 2- Basic graphs to obtain the recursion relations for an Ising model 
on a diamond hierarchical lattice with q = 2 and exchange interactions given 
by the generalized Fibonacci rules, A — > AB, and B — > AA. 

Fig. 3- Second iterates and fixed points (black circles) of the recursion 
relations in the parameter space, xa = tanhi^ versus x B = tanhK B , for 
an Ising model on a diamond hierarchical lattice with q = 2 and exchange 
interactions given by the generalized Fibonacci rules, A — > AB, and B — > 
(irrelevant fluctuations). We draw the stable trivial fixed points (at zero and 
infinite temperatures) and the physical saddle point. The dashed lines and 
the black arrows indicate the flows of the second iterates of the map. The 
uniform model is recovered along the diagonal, xa = %b- The light arrows 
indicate the stable direction in the neighborhood of the physical fixed point. 

Fig. 4- Second iterates and fixed points for the Ising chain: (a) exchange 
interactions according to the rules A — > AB, and B — > AA (irrelevant fluc- 
tuations); (b) exchange interactions according to the rules A — > ABB, and 
B — > AAA (relevant fluctuations). The diagonal solid line corresponds to 
the uniform model. 

Fig. 5- Second iterates and fixed points (black circles) of the recursion 
relations in the parameter space, xa = tanni^ versus x B = tannics, for 
an Ising model on a diamond hierarchical lattice with q = 2 and exchange 
interactions given by the generalized Fibonacci rules, A — > ABB, and B — > 
AAA (relevant fluctuations). We draw the stable trivial fixed points (at zero 
and infinite temperatures) and the unstable node. The uniform model is 
recovered along the diagonal, xa = %b- The dashed lines are indications of 
the flow associated with the second iterates of the map. The light arrows 
indicate the directions of the eigenvectors in the neighborhood of the physical 
fixed point. 
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